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The situation in which the internal conductor of a coaxial transmission line projects through 
a ground-plane without any change in diameter is analyzed with reference to the infinite antenna 
and the parallel-plate antenna. An experimentally meaningful value of the admittance which takes 
into account the geometry of the antenna-to-coaxial line junction is obtained. This is achieved by 
utilizing a mathematical model in which an annular aperture or “gap” is represented by a frill of 
of magnetic current. The field distribution in the aperture is assumed to be of the same form as the 
TEM mode alone in the coaxial line. It is shown that the effect on the admittance of neglecting the 
TM modes in the line at the aperture, is very small. In the case of the parallel-plate antenna, the 
theoretical results are shown to agree remarkably well with experiment. It is also shown how the 
theory may be applied to the problem of finite length cylindrical antennas. 


1. Introduction 


The outstanding problem of cylindrical-antenna theory today is to determine a realistic 
value for the antenna admittance which takes into account the geometry of the feed system. 
In the past, attention has been focused primarily on the evaluation of the current distribution 
on the antenna proper. The nature of the feed system has in general been ignored, and its effect 
obliterated by the use of an infinitesimal “gap.” As is well known, this type of source leads to 
an infinite value for the input susceptance of the antenna when the admittance of a symmetrical 
antenna, say, is defined by the equation, Y=1(0)/V. Attempts to overcome this obviously non- 
physical result have often been directed towards subtracting the singularity in some meaningful 
way (e.g., Chen and Keller, 1962; Duncan, 1962; Rao, 1965). The result of some other investigations 
has been to “‘smooth out“ the singularity (e.g., Duncan and Hinchey, 1960). Some authors have 
bypassed it by using a variational technique (e.g., Storer, 1950; Rao, 1962), and some have defined 
the admittance by the equation, Y=/(d)/V, where d # 0, so that Y is no longer singular (e.g., Infeld, 
1947). These methods are unsatisfactory in that reasonable correlation of theory with experiment 
is either fortuitous, since the results are invariant with respect to the actual aperture size, em- 
pirical, or obtained by way of arbitrary approximations. 

In order to develop a rational and practical theory it is first necessary to decide upon an 
antenna-feed configuration which gives rise to an analytically convenient mathematical model 
without sacrificing practical usefulness. It is also necessary to select the most convenient and 
most widely applicable mathematical method in order that a range of problems encountered in 
practice might be dealt with in a unified way. Finally, it is advantageous if experimental results are 
available for an antenna configuration for which the mathematical analysis is particularly simple in 
order to verify the theory. All these requirements may be met in the following way: First, we use 
the feed geometry of figure 1, the mathematical model of which is shown in figure 2. This model 
was introduced by Albert and Synge (1948). Second, we use the Fourier-transform method. Third, 
in order to verify the theory, we analyze the mathematically simple parallel-plate antenna. 

With regard to the use of the frill model, it is necessary to specify. the tangential electric 
field in the aperture or “gap”’. Thus we neglect all the TM modes in the coaxial line at the aperture 
and take the field distribution to be that of the remaining TEM mode. A plausible argument is 
advanced to show that this approximation should have only a small effect on the resultant value 
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Figure 2. Mathematical model of the antenna. 


FicurE 1. Coaxial line to antenna geometry. 


for the admittance. This argument is reinforced by an analysis of the effect of individual TM 
modes on the admittance. 


Concerning applications, the analysis of the physically unavailable infinite antenna is of 
fundamental significance when the Fourier<transform method is used, for the results therefrom 
may be applied in the solution of problems concerning both parallel-plate and finite-length antennas. 
Thus the first problem considered in this paper is that of the infinite antenna. 


2. The Infinite Antenna with Coaxial Feed 


The antenna under consideration may be taken to be an infinite, perfectly conducting rod 
of radius a mounted vertically above a highly conducting ground plane and fed coaxially from 
below, as shown in figure 1. The cylindrical coordinates ‘(r,6,z) relative to this antenna are also 
shown. Time variation is assumed to be harmonic and proportional to e 


Looking at the problem as a whole shows that there is an essential dichotomy in the type of 
information we want to extract. First, in order to evaluate the radiation characteristics of the sys- 
tem, we must determine the current distribution on the cylindrical part of the antenna as well as 
the total radiation effect of the transmission line. Second, in order to facilitate design of associated 
electronic equipment, we need to know the admittance the antenna presents to the generator. The 
first problem may be called the “exterior” problem and the latter, the “interior” problem. 

This dichotomy suggests that although the two problems are uniquely related, a different 
approach may be warranted in each case. Thus we consider first the interior problem. 
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2.1. The Interior Problem 


Let us effect the division between the interior and exterior problems by means of a fictitious 
surface across the aperture at z=0. Suppose, for the moment, that E, H are known on this surface. 
Then the interior problem is seen to be one of determining the proper combination of incident and 
reflected waves which will give rise to the given E, H on the terminating surface. 

If the usual assumption of perfect conductivity of the coaxial line is made, and if it is also 
assumed that all the fields are axially symmetric and that the only nonzero components of the field 
are E,, Ez, Ho, then He satisfies the wave equation, 

{2 O28 


slat Pere eae eee 2 = 
nor or aoe \He 0, () 


where k?= @pWo€o. (2) 
Once the solution of (1) is known, the remaining field components may be determined from Max- 


well’s equations. Thus, when this is done, and the boundary conditions Ea, z)=0, E.(b, z)=0 
applied, we find that 


E,(r, 2)=— a {Avett— Boe} — : : s Bre*n?Ci(qnr), (3) 
@€or ime) 
Ho(r, j=7 {Avett + Bae + oy Bren?C (gar), (4) 
n=1 
where Cn Gat) = J m( nt) Yo( Gna) — Yn(qnr)Jo(qna). (5) 
The gn are roots of the equation 
Co(qnb) = 0, (6) 
Jn and Y, are Bessel functions, and 
Va qi — F. (7) 
At z=0, we have 
Zo ] =. 
E,(r,0)=— = (Ao— Bo) -——— SS) ynBnCi(Qnr), (8) 
r LWEo awi 
4 s 
Har, )="(to+Bo)+ S BaCilaar) (9) 
n=1 


If we integrate both sides of both these equations with respect to r from r=a, to r=), it is easily 
verified that 


” E,(r,0)dr-=— Ze(4o— Byintb/a), (10) 
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i Hor, 0)dr= (Aot Bo)ln(b/a). (11) 


Therefore, since the admittance of the terminated line is given by 


y= 20 Ee 4] 
Zoln(b/a) Bo + Ag : 


b 
Qn i Hr, 0)dr 


y= —— 2 
Zoln(b/a) | : E,(r, 0)dr oe 


we find that 


This is a most important result for it shows that the admittance is related to the integrals of the tan- 
gential fields in the aperture. 

Let us now assume that the contribution to the admittance from the TM modes in the coaxial 
line at the aperture is small. Thus the tangential electric field in the aperture may be taken to be 


V 


Er, 0)= Tin @Jay’ 


(13) 
where V is the applied voltage across the line. 


A plausible argument in support of this approximation is as follows: First, as r—>a, the E, 
field becomes proportional to 1/r provided ka is not too large. Second, as r— B, it is evident that 
the E, field in the aperture must deviate from the assumed distribution. Indeed there must be an 
edge singularity here as r—> b. However, it is clear from (10) that only the TEM mode in the electric 
field contributes to the integral of FE, in (12). Moreover, it is clear from (8) and (9) that H¢ is “smooth” 
compared with E, because of the factor y, in (8) and its absence in (9). Because of this, and because 
Y is related to the integral of Hg in the aperture so that oscillatory disturbances tend to integrate 
to zero, perturbations from the assumed distribution will give rise to smaller perturbations in the 
admittance. In other words, the main contribution to the result for the admittance is from the TEM 
mode in the coaxial line. The effect of the TM modes should be slight. 

Thus, assuming this argument to be valid (further consideration of the point is deferred until 
section 4) and using (13), we find that the admittance is given by 


2a b 
= Th bla eal Ho(r, 0)dr. (14) 


2.2. The Exterior Problem 


Subject to (13), the exterior problem is completely specified by the following boundary 
conditions: 


E.(a, 2) =0, 4 > 0, 


r ln (b/a) 


Er, = |For a<r<b, 
0,r>6, 
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as well as the usual radiation conditions at infinity. It follows that the only nonzero components of 
the resultant field are E,, Ez, He. 

Evidently we may have recourse to the method of images to simplify the problem. This gives 
rise to the situation shown in figure 2 where the infinite antenna is fed from a frill of magnetic 


current, 
L=— u, xX {E(r, 0+0)—E(r, 0—0)}- (15) 
2V 
. . = | — < 
This gives f rin (Bla)” a<r<b, (16) 
0,r>o6. 


2.3. The Current on the Infinite Antenna 


We are now in a position to determine the current distribution on this antenna. To achieve this 
we observe that the electromagnetic field is determined if the electric and magnetic vector poten- 


tials A and M are known. Thus 


E=-—{VV- +}A~ XM. (17) 


1@[Lo€o 


Because of the symmetry of this problem, it is evident that the only nonzero components of A and 
M are A, and Mg. Consequently 


Baz (Gat ie), fae (rMe), (18) 
_ and A, and Mg are to be determined from 
{2 oF tila =~ poK d2)5(r ~ a) 
rar’ or az ‘ i (19) 


(2 pee ea i|Me =— eb F)(2), (20) 


ror or Pr +32 


where K; is the surface current density on the cylinder. Let us now take the Fourier transform of 
(18), (19), and (20). This transform, where it exists, is defined by the relation, 


¥(r, a) = i _ E(r, 2e- dz. (21) 


Thus we find, for instance, that 


st dr, o) =" (ophuiga)Kolaa). (22) 


Here #(a) is the transform of the total current at any cross section z=constant. Also Ip and Ko 
are modified Bessel functions, g=(a?—k?)"?, and the branch of the square root which is positive 
for large positive a is taken. Moreover, the wavenumber k is assumed to have a small negative 
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imaginary part so that the Fourier transforms exist in the ordinary sense. 
. Solution of the transform of (20) is facilitated by solving for the case of the delta function right- 
hand side with 


Ldr)= 8(r—r’). (23) 
This gives rise to f(r, r’, a), say, and then 
b 

Mr, = i Mr, r', @Ler’)dr’. (24) 


Thus carrying out the indicated calculations and applying the boundary condition that E.(a, z)=0 
in (18), we obtain the result, 


I$ (q) = TRV [x — Kota]. 


Zo In (blag? Ko(qa) = 
Therefore the current on the infinite antenna is given by 
QikV eiaz peice) | 
I2=~— —e | da, 
OF in(la) Ir, | Ko(qa) oe 


where I; is the contour shown in figure 3. Note that this result may be inaccurate for small z 
because, although the effect of the TM modes on the admittance appears to be small, this is not 
necessarily the case with regard to the current. 


FicureE 3. Contour of integration for Fourier inverse. 
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2.4. Numerical Evaluation of the Current 


Numerical evaluation of I(z) from (25) is facilitated by changing the contour of integration to 
T, as shown in figure 4, and making the change of variable, a=— ik Vv?—1. The result is 


I(z)= 


4ikV [= eV 1, ee iy (27) 


~Zoln(bla)Jo yVp—] J2(vka) + Y2(vka) 


We have carried out the above integration using a 1620 digital computer. The result for 
V=%2 V, ka=0.02, and b/a=2.75 is shown in figure 5. Note that we have plotted the amplitude 
term, (kz), in the equation 


I(kz)= Wlkzje~"*, z > 0. (28) 
The corresponding amplitude term in the case of the infinitesimal gap is also shown in this figure. 
Note that the imaginary part of (kz) for the frill source is finite at z=0, whereas that for the infini- 


tesimal gap is logarithmically singular. 


2.5. Calculation of the Admittance 


> 
From (14) we see that the quantity which must be determined is | H,(r, 0)dr. To do this we 
observe that the magnetic field may be determined from . 


w 


Ref¥1, frill model and .5-gop model 


Re[y], A ee 


~ 


Figure 5. Current amplitude on infinite antennas, 
ka= 0.0664, b/a=2.25. 


FicuRE 4. Modified contour of integration. 
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_ 1 4 pe 2s 
icostoes {vv +k fut Vx A. (29) 
For this particular antenna we find that 
1 0A, 
Ho=—ieMo——- (30) 


Therefore in terms of the Fourier transforms of the above quantities, 
f? . f? I 
|. Hae, edr=—io [ aa, adr — 2 {abel )— la, 2)}- (31) 


Using (24), it is easily verified that 


2ikV 


Ks(qb) 
J ete, endr—= ze | in Bla) + FA) (r(gayKe(ad)—Lolab\Kolaa)} | a 


Thus the admittance of the infinite antenna with frill source is given by Y;, where 


vy Dik Ko(qb) _ do, 
ee ERE | . [in (6la)+ Foe IelaadKolab) To(ab)Ko(aa)} | - (33) 


Finally, in a similar manner to the development of (27) for the current, we find that a convenient 
representation for Y; is 


{Jo(vka)Yo(vkb) —Jo(vkb)Yo(vka) }® 
Y= — ik __ [* iJo(vka)Yo(vkb) — Jo(vkb)Yo(vka) }? 1, 
* Zo in a} ie vVv?—1 {Jo(vka) + Y3{vka)} “ C4 


3. The Parallel-Plate Antenna 


Once the Fourier transform of the current on the infinite antenna has been obtained, it is a 
very simple matter to determine the current on a parallel-plate antenna of the configuration shown 
in figure 6. This may be done by constructing a periodic structure of frills of magnetic current excit- 
ing an infinite perfectly conducting rod. If these frills are located atz=+2nh,n=0,1,2 .. ., then 
on the planes z=-(2n— 1h, we get E-=0. The parallel-plate antenna may be identified with one 
cell of this structure, namely, the region 0 < z < A, because the boundary conditions which charac- 
terize the antenna are satisfied. 

Clearly, the current on the periodic structure is given by 


[(2=Ih+2 > In cos (maz/h). (35) 
m1 
But the Im are known to be periodic samples of .4(a), as given by (25), according to the equation, 


1 
In= oh J (mah). (36) 


Thus the current on the parallel-plate antenna is 
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T(z)= 


2riV {Hatkh)_ i}+ 25 coe |, (37) 


Zokh ln (b/a) |. | Holka) Ko(qmka) 


m=1 


mr\? 
where Qn =| (Fr) 4° (38) 


The admittance can be determined from (14) and (32). The result is 


7 oe 
Ho(ka 


ee [ In(b/a) tS | Jk Yate) — J oka) k)| 


Zokh { ln(b/a)}? 


SoS Ge {intbja) + FAG ttaamka\Koldmkb) —Loamk®\Keanka)}}|- 39) 


3.1. Numerical Results 


Equations (37) and (39) were evaluated using a 1620 digital computer. Both series were trun- 
cated at m=20. Values of the current and admittance were calculated for ka =0.0664, b/a=2.25, 
in order that a comparison could be made with some experimental and theoretical results obtained 
by Rao (1962, 1965). 

The current distribution for kh 7/2 is shown in figure 7. Clearly, our theory is satisfactory 
except perhaps near z~0, where the effect of neglecting the TM modes seems to be evident. 
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FicuRE 7, Current on parallel-plate antenna. 
Experimental values (Rama Rao) —§-_ => 
Theoretical values (Otto) ___ x ____. 


Figure 6. Parallel-plate configuration. 
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In table (1), values of both the conductance, Go, and the susceptance, Bo, obtained by the use 
of.(39), are shown for a range of values of kh. Also shown in this table for the purpose of comparison, 
are the corresponding experimental values for the susceptance, Be (Rao, 1965) and also some 
theoretical values. Of these, B, denotes the values due to Rao (1965) when the singularity due to 
a delta-function source is subtracted. The term B, denotes the values obtained when the singularity 
due to the same type of source is avoided by using a variational method (Rao, 1962, 1965). 


TABLE 1. Admittance in mmhos of the parallel-plate antenna, 
(Y =G+iB). 
ka= 0.0664, b/a = 2.25. 


The improvement resulting from the use of our theory is evident from a comparison of these 
results. It is to be noted that Rao’s results, denoted by By, are the same whatever value of b/a is 
taken. Thus the apparent improvement in his results over ours in this case is fortuitous, 


4. The Effect of the TM Modes on the Admittance 


In order to estimate the effect of the TM modes on the admittance, consider the distribution, 


Ci(qmr) 
aCi(qma) In (b/a)’ 


Ex,(7, 0) = (40) 


in which it is evident that 
E,, (a, 0) =E,(a, 0), 


1 


where E roll > » =v in (Bla) ; 


so that the TM modes are roughly normalized with respect to the TEM mode. With this TM-mode 
excitation, the magnetic current on the frill is 


Lg ()=— 2E,_(r, 0). (41) 
Thus the method of sections 2.3 and 2.4 may be used to obtain the result, 


° Dik Ci(qmb)Ko(qb) 1 _ - 
I. HAs adr=F In (bla) Ci(qna)Ko(qa) ate @ x {Io(qa)Ko(qb) I(qb)Ko(qa)} (42) 
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Now let the distribution in the aperture be 
| E(r, 0)=E,,(r, 0) +E,,,(r, 9). (43) 
With this distribution, we can write the result for the admittance in the form, 
Y=Yo+Yn, | (44) 


where Yo is the admittance as calculated using the simple theory and may be identified with the 
first term on the right of (43). The value Ym is the contribution to the admittance from E,,(r, 0). 
Thus, using the technique of section 3, the admittance of the parallel-plate antenna may be deter- 
mined using (42). The result for Ym is, 


Ym =Z.kh{ In (bla)}? C:(pmka) | 2 p31 H@(ka) Fa X {Jo(kb)Yo(ka) — Jo(ka)Yo(kb) } 


27i(b/a) mea 1 H?(kb) 
1 Ko(gnkb) 


+2 
n=1 pit 2 Ko(qnka) 


{Tolqnka)Ko(qnkb) — Lo(qnkb Kolanko) } |, (45) 


where pm=kqm and the qn are defined by (6). 


The term Y, was evaluated for ka= 0.0664, b/a=2.25, and m=1, 2, and 3, for a few values 
of kh. The results are shown in table (2) along with Yo for comparison. Also shown is a series of 
values of Go, Bo, Gi, B:, for some of the values of kA given in table (1). Clearly, the contribution 
from the Yn is small, thus verifying the argument put forward in section 2.1. Apparently this result 
is due to the factor 1/(p2+ gq?) in (45), because pm ~ ma/{ka(b/a — 1)}. Thus for small ka, Ym will 


be small. Also note that the purturbation to the aperture distribution implied in (43) is quite vio- 
lent. Physically, one would expect the actual deviation to be much smaller. Thus the resultant 
effect of the TM modes is probably much smaller than that indicated by table (2). This raises the 
question of the origin of the significant constant difference of about 1 mmho between Bo and B, 
in table (1). The point appears to be worthy of further consideration. 


TABLE 2. Effect of TM modes on the admittance of the parallel- 
plate antenna, (Y =G + iB). 
All admittance values in mmhos, ka = 0.0664, 5/a= 2.25. 


1 
2 
3 
1 
2 
3 
1 
2 
3 
1 
1 
1 
1 
1 
1 
1 
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5. The Admittance of a Finite Antenna 


Einarsson’s (1966) method may be used to show that the current on the cylindrical surface of 
a finite antenna of length h, driven by a frill source, is given by 


I2)=f2)— i ee'{Ielz— h—z’)+L(zth+z)}de', (46) 


where Jz) is the current on an infinite antenna with frill source, J-(z) is the current on a similar 
antenna with delta-function source, and e(z) is the z component of the electric field on the cylinders 
r=a,z>h. 

An approximate value for the admittance may be found if we observe that the right of (46) 
consists of an outgoing current wave and its “reflections” from the ends of the antenna. The 
Hg field due to these reflections is approximately proportional to 1/r over the aperture provided 
the size of the aperture is small compared with the length of the antenna. Thus the admittance of 
the finite antenna is 


Y=Y,-2 | e(z'We(h + 2"\d2’ , (47) 


where Y; is the admittance of the infinite antenna with frill source as given by (34). 


6. Conclusion 


We have shown that the use of a frill model for an antenna driven from a coaxial line, along 
with the approximation that the tangential electric field in the aperture is of the same form as the 
TEM mode in the line, leads to a theory which is easy to apply and gives results which correlate 
very well with experiment. In the case of the determination of the admittance, the extent of the 
above approximation has been assessed by evaluating the contribution to the admittance of individ- 
ual TM modes in the aperture. This contribution was shown to be very small for small ka. Thus 
the theory has the distinct advantage over those mentioned in the introduction, that it gives a 
rational and accurate way of taking into account the dimensions of the coaxial aperture. 
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